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A strongly coupled quark-gluon plasma (QGP) of heavy constituent quasiparticles is studied by 
a path-integral Monte-Carlo method, which improves the corresponding classical simulations by 
extending them to the quantum regime. It is shown that this method is able to reproduce the 
lattice equation of state and also yields valuable insight into the internal structure of the QGP. The 
results indicate that the QGP reveals liquid-like rather than gas-like properties. At temperatures 
just above the critical one it was found that bound quark-antiquark states still survive. These states 
are bound by effective string-like forces. Quantum effects turned out to be of prime importance in 
these simulations. 



I. INTRODUCTION 



Investigation of properties of the quark-gluon plasma (QGP) is one of the main challenges of strong-interaction 
physics, both theoretically and experimentally. Many features of this matter were experimentally discovered at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven. The most striking result, obtained from analysis of these 
experimental data [H, is that the deconfined quark-gluon matter behaves as an almost perfect fluid rather than a 
perfect g it could be expected from the asymptotic freedom. 

There are various approaches to studying QGP. Each approach has its advantages and disadvantages. The most 
fundamental way to compute properties of the strongly interacting matter is provided by the lattice QCD [3-13]. In- 
terpretation of these very complicated computations requires application of various QCD motivated, albeit schematic, 
models simulating various aspects of the full theory. Moreover, such models are needed in cases when the lattice QCD 
fails, e.g. at large baryon chemical potentials and out of equilibrium. 

A semi-classical approximation, based on a point like quasi-particle picture has been introduced in [5]. It is 
expected that the main features of non-Abelian plasmas can be understood in simple semi-classical terms without the 
difficulties inherent to a full quantum field theoretical analysis. Independently the same ideas were implemented in 
terms of molecular dynamics (MD) [6j . Recently this MD approach was further developed in a series of works 0, ll| . 
The MD allowed one to treat soft processes in the QGP which are not accessible by perturbative means. 

A strongly correlated behavior of the QGP is expected to show up in long-ranged spatial correlations of quarks 
and gluons which, in fact, may give rise to liquid-like and, possibly, solid-like structures. This expectation is based 
on a very similar behavior observed in electrodynamic plasmas f?', 19^ . This similarity was exploited to formulate a 
classical non-relativistic model of a color Coulomb interacting QGP [7,] which was numerically analyzed by classical 
MD simulations. Quantum effects were either neglected or included phenomenologically via a short-range repulsive 
correction to the pair potential. Such a rough model may become a critical issue at high densities. For temperatures 
and densities of the QGP considered in Ref. Q these effects are very important as the quasiparticle thermal wave 
length is of order the average interparticle distance. 

In this contribution we extend previous classical nonrelativistic simulations 0] based on a color Coulomb interac- 
tion to the quantum regime. We develop an approach based on path integral Monte Carlo (PIMC) simulations of 
the strongly coupled QGP which self-consistently takes into account the Fermi (Bose) statistics of quarks (gluons). 
Following an idea of Kelbg [10.] , quantum corrections to the pair potential are rigorously derived [llj . This method has 
been successfully applied to strongly coupled electrodynamic plasmas [H, [isl . Examples are partially-ionized dense 
hydrogen plasmas, where liquid- like and crystalline behavior was observed [l4l [isj . Moreover, also partial ionization 
effects and pressure ionization could be studied from first principles p^ . The same methods have been also applied 
to electron-hole plasmas in semiconductors [iTlllsj. including excitonic bound states, which have many similarities to 
the QGP due to smaller mass differences as compared to electron-ion plasmas. 

First results of applications of the PIMC method to the nonideal QGP have already been briefly reported in [19,] . 
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II. BASICS OF THE MODEL 

Our model is based on a resummation technique and lattice simulations for dressed quarks, antiquarks and giuons 
interacting via the color Coulomb potential. The assumptions of the model are similar to those of Ref. Q : 

I. All color quasi-particles are heavy, i.e. their mass (m) is higher than the mean kinetic energy per particle. For 

instance, at zero net-baryon density it amounts to to > T, where T is a temperature. Therefore these particles 
move non-relativistically. This assumption is based on the analysis of lattice data [1^ [2l| . 

II. Since the particles are non-relativistic, interparticle interaction is dominated by a color-electric Coulomb potential, 

see Eq. ([1]). Magnetic effects are neglected as sub- leading ones. 

III. Relying on the fact that the color representations are large, the color operators are substituted by their average 

values, i.e. by classical color vectors, the time evolution of which is described by Wong's dynamics [23 |. 

The quality of these approximations and their limitations were discussed in Ref. Thus, this model requires the 
following quantities as an input: 

1. the quasiparticle mass, to, and 

2. the coupling constant . 

All the input quantities should be deduced from the lattice data or from an appropriate model simulating these data. 

Thus, we consider a three-component QGP consisting of a number of dressed quarks (iVg), antiquarks {Nq) and 
giuons {Ng) represented by quasi-particles. In thermal equilibrium the average values of these numbers can be found in 
the grand canonical ensemble defined by the temperature-dependent Hamiltonian, which can be written as H = K+U. 
The kinetic and color Coulomb interaction energy of the quasi-particles are 

2 2^ An\r,-rA ' ^ ' 



2TOi(T, fj.q) 



2J 



Here the Qi denote the Wong's color variable (8-vector in the SU{3) group), T is the temperature and /Zg is the quark 
chemical potential. In fact, the quasi-particle mass and the coupling constant, as deduced from the lattice data, are 
functions of T and, in general, ^q. Moreover, g^ is a function of distance r, which produces a linearly rising potential 
at large r [23;]. 

The thermodynamic properties in the grand canonical ensemble with given temperature T, chemical potential /ig 
and fixed volume V are fully described by the grand partition function 

Z{^^q,p,V)^ Y: ^"p(M^-^;?)/^) ^ / drdQp{r,Q,a;Nq,Nq,N,-J), (2) 

where p{r,Q,a; Nq, Ng, Ng] l3) denotes the diagonal matrix elements of the density operator p = exp(— and 
/3 = 1/T. Here ct, r and Q denote the spin, spatial and color degrees of freedom of all quarks, antiquarks and giuons 
in the ensemble, respectively. Correspondingly, the a summation and integration drdQ run over all individual degrees 
of freedom of the particles. Since the masses and the coupling constant depend on the temperature and chemical 
potential, special care should be taken to preserve thermodynamical consistency of this approach. In order to preserve 
the thermodynamical consistency, thermodynamic functions such as pressure, entropy, S, baryon number, Nb, 
and internal energy, E, should be calculated through respective derivatives of the logarithm of the partition function 

P ^ d{TliiZ)/dV, S^d(TlnZ)/dT, iVs = (l/3)9(rinZ)/9^g, E ^ -PV + TS + SfigNe- (3) 

This is a conventional way of maintaining the thermodynamical consistency in approaches of the Ginzburg-Landau 
type as they are applied in high-energy physics. 

The exact density matrix of interacting quantum systems can be constructed using a path integral approach [23 . [25j 
based on the operator identity e~^^ = e~^^^ ■ e'^^^ . . . e~^^^ , where the r.h.s. contains n + 1 identical factors 
with A/3 = (i/in + 1). which allows us to rewrite^ the integral in Eq. ([2]) 

^ f dr'^'^UQ^''^ p{q(''\Q(''\<T-Nq,Nq,Ng;!3)= f drWdQ(o)dr«dQ(i) . . . dr(")dQ(") p« • . . . 



^ For the sake of notation convenience, we ascribe superscript to the original variables. 
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p p- p 

JTq JTq JTg 

dQ^°^ dr(°) dr(i) . . . dr p(r , r , . . . r ; g(0) ; iV, , iV,-, iVg ; /3) . (4) 



The spin gives rise to the spin part of the density matrix (S) with exchange effects accounted for by the permu- 
tation operators Pq, Pq and Pg acting on the quark, antiquark and ghion spatial r("+^) and color coordi- 
nates, as well as on the spin projections a' . The sum runs over all permutations with parity Kp^ and up-. In 
Eq. (IH) the index / = + 1 labels the off-diagonal density matrices p"^'^ = p (r^'"^), Qf'"^); r^'), Q^'^; A/3) « 

^j,(i-i)|g-A/3ff|^(0^5^(Q(i-i) _ q(0)^ where 5^{Q'^^~^^ - is a deha-function at e ^ 0. Accordingly each a particle 

is represented by a set of n -I- 1 coordinates ("beads"), i.e. by {n + 1) 3-dimensional vectors {ri""*, . . .ri""*} and a 
8-dimensional color vector Q'^^'^ in the S'C/(3) group, since all beads are characterized by the same color charge. 

The main advantage of decomposition Q is that it allows us to use a semi-classical approximation for density 
matrices p^^\ which is applicable due to smallness of artificially introduced factor l/{n+ 1). This parameter makes 
the thermal wavelength AAq = •^/27rA/3/m^ of a bead of type a {a = q,q,g), smaller then a characteristic scale of 
variation of the potential energy. In the high-temperature limit can be approximated by a product of two-particle 
density matrices. Generalizing the electrodynamic plasma results to the case of an additional bosonic species (i.e. 
gluons), we write 

p(rl»>,r"l,...rWi<3<»'iJV„JV5,«,;/!) 

- E ^^ °;S;f^i:;?»f" ii^-'ik.--"-' iij"'°ii.detii j-°iu. n n c 

s,k ■^l -^9 1=1 p=l 

where N = Nq + Nq + Ng, s and k are numbers of quarks and antiquarks, respectively, with the same spin projection, 
Xa = \/2Tr/3 /ma is a thermal wavelength of an a particle, = Nal/[s\{Na — s)!], the antisymmetrization and sym- 
metrization are taken into account by the symbols "det" and "per" denoting the determinant and permanent, respec- 



tively. Functions (f>pp = exp 



(0 ^ 



and matrix elements = exp 



(rf-r(°))+y^ /AXi\6,{Q, 



Qo), where t and o are particle's indexes, are expressed in terms of distances (yi^\ . . . , yi"-*) and dimensionless 
distances (^i^"*, . . . , ^i""*) between neighboring beads of an a particle, defined as = r'^^ + ya\ {I > 0), and 
Ua^ = AXaJ2k=i^a'^ ■ Noticc that the ||0"'°|| matrix consists of three nonzero blocks related to quarks, 
antiquarks, and gluons, 1 |giue- The density matrix (O has been transformed to a form which does not 

contain an explicit sum over permutations. Let us stress that the determinants depend also on the color variables. 
In Eq. ([5]) the total color interaction energy 

, n+l 

p^t 1=1 

is defined in terms of off-diagonal two-particle effective quantum potential which is obtained by expanding the 
two-particle density matrix ppt up to the first order in small parameter l/(n + 1): 

Pptir,r',Qp,Q„Al3)^p%{r,r',Qp^Q,^A(3)~ tdr [ ^^^^ ^^dHW'lT, ^,q){Qp\Q,) 



4^|r"|AA2,V^(l-T) 

{-^^;=^)^^v{-^^^ (7) 



where r — Vp rt, r' = r'^ — r[, AXpt = ^27rA/3/mpf , nipt = fnpmt/{nip + toj), and is the two-particle density 
matrix of the ideal gas. The result for the diagonal color Kelbg potential can be written as 

$^*(r,r,Qp,O0«^^%g^^^^ {l-e-?*+V^:.p*[l-erf(xpO]}, (8) 



where Xpt = \rp — rt\/ AXpt- Here the function g^{T,pq) = g'^{r" ,T.p,q), resulting from averaging of the initial 
g^{r" , T, Pq) over relevant distances of order A Apt, plays the role of an effective coupling constant. Note that the color 
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Kelbg potential approaches the color Coulomb potential at distances larger than A Apt. What is of prime importance, 
the color Kelbg potential is finite at zero distance, thus removing in a natural way the classical divergences and making 
any artificial cut-offs obsolete. This potential is straightforward generalizations of the corresponding potentials of 
electrodynamic plasmas [1^ . The off-diagonal elements of the effective interaction are approximated by the diagonal 
one by means of $P*(r, r'; , Qp, Qt) « [$f*(r, r, Qp, Qt) + $P*(r', r', Qp, Qt)]/2. 

The described path- integral representation of the density matrix is exact in the limit n — > oo. For any finite number 
n, the error of the above approximations for the whole product on the r.h.s. of Eq. ^ is of the order l/{n + 1) 
whereas the error of each is of the order -I- 1)^, as it was shown in Ref. [l3| . 

The main contribution to the partition function comes from configurations in which the 'size' of the cloud of beads 
of quasiparticles is of the order of their thermal wavelength Aa whereas characteristic distances between beads of each 
quasiparticle are of the order of AA^ . 



III. TESTING THE METHOD WITHIN THE CANONICAL ENSEMBLE 

To test the developed approach we consider the QGP only at zero baryon density and further simplify the model 
by additional approximations, similarly to Ref. 01 : 

IV We replace the grand canonical ensemble by a canonical one. The thermodynamic properties in the canonical 

ensemble with given temperature T and fixed volume V are fully described by the density operator p = e^^^ 
with the partition function defined as follows 

Z{N,,N,,Ng,V;f3) - ^ , E / drdQ p{r,Q,a; f3), (9) 

« 9 9 o- ^ 

with Nq = Ng and hence Nb — 0. In order to preserve the thermodynamical consistency of this formulation, 
thermodynamic quantities should be calculated through respective derivatives of the logarithm of the partition 
function similarly to that in Eq. ([3]) with the exception that now Na are indepenent variables. 

V Since the masses of quarks of different fiavors extracted from lattice data are very similar, we do not distinguish 

between quark flavors. Moreover, we take the quark and gluon quasi-particle masses being equal because their 
values deduced from the lattice data [1^ [2l| are very close. 

VI Because of the equality of masses and approximate equality of number of degrees of freedom of quarks, antiquarks 

and gluons, we assume that these species are equally represented in the system: Nq = Nq = Ng. 

VII For the sake of technical simplicity, the SU(3) color group is replaced by SU(2). 
Thus, this simplified model requires an additional quantity as an input: 

3. the density of quasi-particles {Nq + Nq + Ng)/V ~ n{T) as a function of the temperature. 

Although this density is unknown from the QCD lattice calculations and we use it as a fit parameter, it is very 
important to partially overcome constrains of the above simplifications. First, it concerns the use of the SU(2) color 
group, which first of all reduces the degeneracy factors of the quark and gluon states, as compared to the SU(3) case, 
and thereby reduces pressure and all other thermodynamic quantities. A proper fit of the density allows us to remedy 
this deficiency of the normalization. Second, in fact we consider the system of single quark flavor, i.e. all quarks are 
identical, which also reduces the normalization of all thermodynamic quantities. The density fit cures the deficiency 
of this normalization, though the excessive anticorrelation of quarks remains. 

Ideally the parameters of the model should be deduced from the QCD lattice data. However, presently this task is 
still quite ambiguous. Therefore, in the present simulations we take a possible (maybe, not the most reliable) set of 
parameters. Following Refs. 0, HH, the parametrization of the quasi-particle mass is taken in the form 

m{T)/Tc = 0.9/{T/Tc - 1) + 3.45 -I- OAT/Tc (10) 

where Tc = 175 MeV is the critical temperature. This parametrization fits the quark mass at two values of temperature 
obtained in the lattice calculations [20j. According to [l^l the masses ar e q uite large: iriq/T ~ 4 and TUg/T ~ 3.5. 
These are essentially larger than masses required for quasi-particle fits [13, of lattice thermodynamic properties 
of the QGP: mq/T ~ 1 ^ 2 and rUg/T ~ 1.5 3. Moreover, the pole quark mass iriq/T ~ 0.8 was reported in recent 
work |29|], as deduced from lattice calculations. Nevertheless, in spite of the fact that it obviously produces too 
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high masses, we use parametrization in order to be compatible with the input of classical MD of Ref. [7|. The 
T-dependence of this mass is illustrated in Fig. [T] (left panel) . 

The coupling constant used in the simulations is displayed in Fig. [T] as well. From the point of view of the QCD 
phenomenology [3(?| it is too high at low energies, i.e. at T /T^ ~ 1 2. However, the high values of masses require 
such a large value of g^, e.g., to be consistent with the HTL results for the quasi-particle masses. The large value of 

is also required to simulate larger values of Casimirs (defining the normalization of the color vectors) in the 5*^7 (3) 
group as compared to the SU{2) one used here. Moreover, such high 5^ are not inconsistent with the lattice data [31|. 

The density of quasi-particles, which is additionally required within the canonical-ensemble approach, was chosen 
on the condition of the best agreement of the calculated pressure with the corresponding lattice result, see Fig. [T] (right 
panel). It was taken to be n{T) = 0.24T^. From the first glance, it is a very low density. For example, in the classical 
simulations of Ref. Q it was taken as n{T)/T^ = 6.3, which corresponds to the density of an ideal gas of massless 
quarks, antiquarks and gluons. Since the quasi-particles are very heavy in the present model (as well as in that of 
Ref. [3]), the latter density looks unrealistically high. Even in quasi-particle models [13, HI], where the masses are 
lower, the density turns out to be n{T)/T^ « 1.4. Since Eq. ([T0|) gives even larger masses than those in Refs. [13, HI] 
and in view of the adopted large coupling, the chosen value of nlT) does not look too unrealistic. The T-dependence 
of n(r) is presented in Fig. [T] (left panel) in terms of the mean interparticle distance rgiT) in units of cr = l/Tc =1-1 
fm (Wigner-Seitz radius). Notice that the present choice of n{T) corresponds to the relation Trs{T) = 1. 

Thus, although the chosen set of parameters is still debatable, it is somehow self-consistent. In the future we are 
going to get rid of the n(T) parameter, by applying the grand-canonical approach, and by using more moderate (and 
maybe realistic) sets of parameters. 

Calculation of the equation of state (right panel of Fig. [T]) was used to optimize the parameters of the model in order 
to proceed to predictions of other properties concerning the internal structure and in the future also non-equilibrium 
dynamics of the QGP. The plasma coupling parameter, F defined as a ratio of the average potential to average kinetic 
energy, is also presented in Fig. [1] (left panel). It turns out to be of the order of unity which indicates that the QGP 
is a strongly coupled Coulomb liquid rather than a gas. In the studied temperature range, 1 < T/Tc < 3, the QGP 
is, in fact, quantum degenerate, since the degeneracy parameter Xa = "-aA^ varies from 0.1 to 2. 
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FIG. 1: Left panel: Temperature dependence of the model input quantities and the plasma coupling parameter F. Right panel: 
Equation of state (pressure versus temperature) of the QGP from PIMC simulations compared to lattice data of Refs. 0, 13| • 



Let us now consider the spatial arrangement of the quasiparticles in the QGP by studying the pair distribution 
functions (PDF's) gabir)- They give the probability density to find a pair of particles of types a and 6 at a certain 
distance r from each other and are defined as 



gab{Ri — R2 



1 



ZN,\N,\Ng\ ^ 



drdQ S{Ri - rj)5(i?2 - r^)p(r, Q, a; /3). 



(11) 



The PDF's depend only on the difference of cordinates because of the translational invariance of the system. In a 
non-interacting classical system, gabif) = 1, whereas interactions and quantum statistics result in a re-distribution of 
the particles. Results for the PDF's at temperature T/Tc — 3 are shown in the top panels of Fig. [2l 

At large distances, r/(j> 0.5, all PDF's of identical particles (top left panel of Fig. coincide, approaching unity. 
At small distances, the gluon PDF increases monotonically when the distance goes to zero, while the PDF of quarks 
(and antiquarks) exhibits a broad maximum. In the present conditions, the thermal wavelength A approximately 
equals 0.37(j, i.e. the difference starts to appear at distances of the order of A. The enhanced population of low 
distance states of gluons is due to Bose statistics and the color-Coulomb attraction. In contrast, the depletion of 
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FIG. 2: Pair distribution functions (upper panels) and color pair distribution functions (lower panels) of identical (left panels) 
and different (right panels) quasiparticles at temperature T /Tc = o = 1/Tc — 1.1 fm. 

the small distance range for quarks is a consequence of the Pauli principle. In an ideal Fermi gas g{r) equals zero 
for particles with the same spin projections and colors, while for particles with different colors and/or opposite spins 
the PDF equals unity in the limit r 0. As a consequence, the spin and color averaged PDF approaches 0.5. 
Such a low-distance behavior is also observed in a nonideal dense astrophysical electron-ion plasma and in a nonideal 
electron- hole plasmas in semiconductors [13, . The depletion of the probability of quasiparticles at small distances 
results in its enhancement at intermediate distances. This is the reason for the corresponding PDF maxima. 

At small distances, r < 0.3a, a strong increase is observed in all PDF's of particles of different type (top right panel 
of Fig. [2]), which resembles the behavior of the gluon-gluon PDF. This increase is a clear manifestation of an effective 
pair attraction of quarks and antiquarks as well as quarks (antiquarks) and gluons. This attraction suggests that the 
color vectors of nearest neighbors of any type are anti-parallel. If this explanation is correct can be verified by means 
of color pair distribution functions (CPDF) defined as 

CabiRi - R2) = zn/n^IN^I ^ / ^'"^^ (Qi\Q2)mi - r1)S{R2 - rl)p{r, Q, a; /3), (12) 

which are shown in the lower panels of Fig. [2] All CPDF's turn out to be negative at small distances, indicating 
anti-parallel orientation of the color vectors of neighboring quasiparticles. The minimum of Cqq close to r = 0.2t7 
corresponds to the maximum observed in gqq. The deep minimum in the gluon CPDF at small distances results from 
the Bose statistics and complies with the high maximum of the gluon PDF ggg. 

Thus, at T/Tc = 3 we observe signs of a spatial ordering, cf. the peak of the quark PDF around r/tr = 0.1 — 0.2, 
which may be interpreted as emergence of liquid-like behavior of the QGP. The QGP lowers its total energy by 
minimizing the color Coulomb interaction energy via a spontaneous "anti-ferromagnetic" ordering of color vectors. 
This gives rise to a clustering of quarks, antiquarks and gluons. 

Fig. [3] presents PDF's of the identical particles for two temperatures T = l.lTc and T = 2Tc (left panels). 
These PDF's can be formed either by correlated scattering states or by bound states of quasiparticles, depending 
on the relative fractions of these states. In a strict sense, however, there is no clear subdivision into bound and free 
"components" due to the mutual overlap of the quasiparticle clouds. In addition, there exists no rigorous criterion 
for a bound state at high densities due to the strong effect of the surrounding plasma. Nevertheless, a rough estimate 
of the fraction of quasiparticle bound states can be obtained by the following reasonings. The product r^gab{r) has 
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FIG. 3: Pair distribution functions at two different temperatures T (left) and quark-antiquark PDF multiplied by distance 
squared (right). 

the meaning of a probability to find a pair of quasipartices at a distance r from each other. On the other hand, the 
corresponding quantum mechanical probability is the product of and the two-particle Slater sum 

E,, = 8^3/2_^3^ ^ |^^(^)|2 exp(-/3i?„) = S^, + Y^l,, (13) 

Q 

where and (r) are the energy (without center of mass energy) and the wave function of a quasiparticle pair, 

respectively, and \ab — \J i-Kh^ j3{ma + mb) / (mamb) ■ Ea6 is, in essence, the diagonal part of the corresponding density 

matrix. The summation runs over all states a of the discrete (Eff,) and continuous (E^^) spectrum. 

At temperatures smaller than the binding energy and distances smaller than or of the order of several bound 
state radii, the main contribution to the Slater sum comes from bound states. In the electromagnetic plasma it was 
found that the product r^Ef^ is sharply peaked at distances around the Bohr radius in this case. Similarly, at low 
temperature, r^gqqir) forms a pronounced maximum near r = 0.2 fm which can be interpreted as the radius of a 
bound qq pair. Thus, our calculations support the existence of bound states of medium-modified (massive) quarks 
and gluons at moderate temperatures, i.e. just above Tc, proposed in Ref. [s^ and later in Refs. [13,^ based on 
results from lattice QCD calculations of spectral functions |3l, [s^]- With the temperature rise these bound states 
dissolve much faster than it was assumed in [sj, HI], which complies with the analysis of Ref. Indeed, at the 

temperature of T = 2Tc the bound states completely disappear (see right lower panel of Fig. [3]). 

Interesting observations can be done from the analysis of the potential of average force (PAF) defined as the 
logarithm of the related PDF, Uab{r,T) = ~T\ngab{r,T). This definition is motivated by the PDF virial expansion 
in terms of bare potential (like color Kelbg potential). Near the QGP phase transition the PAF is a linear function 
at distances smaller than the bound state radius. This suggests that the bound states are bound by a string-like 
forces. At larger distances the PAF can be very well approximated by an exponentially screened Coulomb potential 
(Yukawa-type potential) like that in the electromagnetic plasma. 

IV. CONCLUSION 

Quantum Monte Carlo simulations based on the quasiparticle picture of the QGP are able to reproduce the lattice 
equation of state (even near the critical temperature) and also yield valuable insight into the internal structure of 
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the QGP. Our results indicate that the QGP reveals liquid-like (rather than gas-like) properties even at the highest 
considered temperature of ST^. At temperatures just above we have found that bound quark- antiquark states still 
survive. These states are bound by effective string-like forces. Quantum effects turned out to be of prime importance 
in these simulations. 

Our analysis is still too simplified and incomplete. It is still confined only to the case of zero baryon chemical 
potential. The input of the model also requires refinement. Work on these problems is in progress. We have also 
performed first simulations of dynamic properties of the QGP based on quantum Wigner dynamics. In particular, 
these allow us to deduce the viscosity of the QGP. However, the brief format of the present contribution does not 
allow us to report on the respective results. 

We acknowledge stimulating discussions with D. Blaschke, M. Bleicher, R. Bock, B. Friman, C. Ewerz, D. Rischke, 
and H. Stoecker. Y.I. was partially supported by the Deutsche Forschungsgemeinschaft (DFG projects 436 RUS 
II3/558/0-3 and WA 43I/8-I), the RFBR grant 09-02-91331 NNIO_a, and grant NS-7235. 2010.2. 
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